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Abstract 

The role of neutrinos in big bang nucleosynthesis is reviewed. Neutrino oscillations 
in the early universe both in resonance and non-resonance case are briefly discussed. 
BBN and supernova limits on heavy sterile neutrinos with the mass in 10-100 MeV 
range are presented and compared with direct experimental bounds. 



1 Basics of BBN 



Big Bang Nucleosynthesis is known as one of the three sohd pillars of the Standard 

Cosmological Model (SCM); the other two are the Cosmic Microwave Background 
Radiation (CMBR) and General Relativity (GR). BBN describes creation of light 
elements in the early hot universe in a very (or rather, depending upon the attitude) 
good agreement with astronomical observations. The theory predicts mass fraction 
of '^He (with respect to total baryonic matter) about 25%, the relative number 
densities of deuterium, ^iJ, and ^He at the level a fewxlO"^ each, and a fewxlO"^'^ 
of '^Li. The predictions span 9 orders of magnitude and well fit the data. 

The characteristic temperature and time scales for the processes of light 
element formations are respectively T = (~ MeV) — 70 keV and t = 0.1 — 10^ sec. 
The relation between the time and temperature (the coohng rate) at this stage is 
given by 

/lO 7^\^^^ 

(i/sec)(T/MeV)2 = 0.74 ( 1 (1) 

Where describes the particle content of the primeval plasma. In the SCM it is 
equal to 10.75 with the following contributions: 2 from photons, 7/2 from electron- 
positron pairs, and 3-7/4 from three neutrino flavors. Any other form of energy 
present in the plasma is parametrized by the contribution into (7* as additional 
neutrino flavors, 5g^ = (7/4) (iV,^ — 3), though these hypothetical forms of energy 
could be either in the form of massive particles, or vacuum-like energy, or anything 
else. Because the cosmological cooling rate depends upon g^, BBN is sensitive to 
any form of energy from MeV down to tens of keV range. 

There are of course some baryons in the plasma, building material for light 
nuclei, but they do not have a noticeable contribution into the total energy density 
because they are quite rare. Their number density, ub, is characterized by the 
parameter 

7710 = — = a few X 10-^*^ (2) 

As we see in what follows, the output of light elements produced at BBN depends 
upon riiQ. It is the only unknown parameter of the standard theory. In fact, the 
data on light element abundances are used to determine riiQ. Till the last year 
it was the only reasonably accurate way. Direct astronomical observations of the 
fraction of baryonic matter in the universe give about 10% of the necessary amount 
of baryons. A dominant part of baryons is not visible. This year the measurements 



of the relative heights of the first and second peaks in the angular fluctuations of 
CMBR permitted to determine rjiQ independently EP : 



ryf/^^« = 5.7±1.0 (3) 

This value is rather close to ?7io determined from BEN (see below). 

Production of light elements crucially depends upon the available number 
of neutrons. The latter is determined by the competition between the rate of the 
reactions of neutron-proton transformation 

n + z/e p + e~, n + e'^ ^ p + u (4) 

and the cosmological expansion rate, 



H = 5.44 J-^ — (5) 
V 10.75 mpi 

where rripi = 1.22 ■ 10^^ GeV is the Planck mass. The rate of the reactions is 
Pnp ~ CjpT^ and at T > 0.7 MeV reactions are faster than expansion and neutron-to- 
proton ratio follows equilibrium value, {n/p)eq = exp(— Am/T), where Gp = 1.166 x 
10^^ GeV~^ is the Fermi coupling constant and Am = 1.293 MeV is the neutron- 
proton mass difference. When the temperature drops below ~ 0.7 MeV the reaction 
become slow in comparison with expansion and the ratio n/p would remain constant 
(frozen) if not the neutron decay with the life-time r„ = 886.7 sec. As a result of 
the decay the neutron-to-proton ratio slowly decreases as exp(— i(:/r„). It goes on till 
the temperature drops down to the temperature of nucleosynthesis, Tbbn = 60 — 70 
keV, when formation of light elements abruptly and quickly begins. Practically all 
neutrons end their lives in "^He since the latter has the largest binding energy. Mass 
fraction of "^He rather weakly, logarithmically, depends upon the number density of 
baryons, rjiQ. This quantity determines the temperature of BEN and respectively 
the number of neutrons that survived to this moment. One can see that Tppj^j is a 
very mild function of rjiQ. On the other hand, helium-3 and especially deuterium are 
very sensitive to ?7io because these elements disappear in two-body baryon collisions 
whose probability is proportional to the number density of baryons. Ey this reason 
the amount of primordial deuterium serves as "baryometer" , measuring rjiQ. The 
calculated abundances of light elements as functions of rjiQ are presented in figs. |1|, 
0, ^. taken from the review ^ . 

According to the recent measurements the observational data on primordial 
mass fraction of ^He is concentrated near two centers Yp = 0.234 — 0.238 @' 0' 
and Yp = 0.244 0). The first value corresponds to r/io ~ 2, while the other is close 
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Figure 1: The predicted ^He abundance (solid curve) and the 2a theoretical uncer- 
tainty @) . The horizontal lines show the range indicated by the observational data. 



to 4-5. The analysis of ref. was based on the data of the work 0) but with an 
account of ionization corrections. 

Recent measurements of the number fraction of deuterium, D/H, in high 
red-shift Lyman-alpha clouds which presumably are not contaminated by stellar pro- 
cesses and have its primordial quantity not diminished by burning is stars, give the 
results between 4 and 1.5 in units of 10~ The corresponding 

values of rjiQ are between 4.5 and 8.5. There are also several measurements of very 
high deuterium values, 15-20 in the same units 0) , for the references to earlier 
papers see e.g. However it was argued recently in ref. 0) that the system at 

red-shift z = 0.7, where the high deuterium value has been observed, has a com- 
plex velocity field and cannot be used for identification of deuterium because the 
deuterium line corresponds to velocity -81 km/sec and can be easily mimicked by 
the turbulent velocity field. It is concluded in the paper that high deuterium is 
practically excluded. 
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Figure 2: The predicted D/H abundance (solid curve) and the 2a theoretical uncer- 
tainty 0) . The horizontal lines show the range indicated by the observational data 
for both the high D/H (upper two lines ) and low D/H (lower two lines). 

2 Role of neutrinos in BBN 

There are several physical effects that make neutrinos important for the nucleosyn- 
thesis. First, as we have already mentioned above, neutrinos contributes to the 
total cosmological energy density and through it to the cooling rate, see eq. (0). 
The freezing temperature of ra/p-transformation depends on g^, as Tfr ~ gl^^ and 
the larger is g^,, the larger becomes ra/p-ratio and more ^He is produced. Another 
effect is a decrease of the time interval which is necessary to reach the nucleosynthe- 
sis temperature Tbbn (the latter does not depend on g^). Both effects work in the 
same direction and result in an increase of the mass fraction of ^He. An addition 
of one extra neutrino species leads to the increase of Yp by approximately 5%. For 
deuterium the effect is of the same sign and stronger, see figs. |^. 

The energy distributions of neutrinos are assumed to have the Fermi-Dirac 
equilibrium form: 




[l + exp(i?/T-e,)] 



(6) 



where E is the neutrino energy, and is dimensionless chemical potential, = 
fij/T, and j denotes neutrino flavour, j = e, fj,, r. In the standard model is assumed 
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Figure 3: The predicted "^Li abundance (solid curve) and the 2cr theoretical uncer- 
tainty . The horizontal lines show the range indicated by the observational data. 

that lepton asymmetry of neutrinos is negligibly small or in other words = 0. 
The impact of chemical potentials on light elements abundances is two- fold. First, 
degenerate neutrinos, i.e. those with a non-zero ^, have a larger energy density, 
corresponding to: 

= 10.75 

Thus a non-zero ^ results in an increase of ^He and deuterium. 

The second effect is operative only for electronic neutrinos because the 
latter are directly involved in the reactions (^. If z/g are degenerate, then the 
equilibrium n/p ratio is shifted by the factor exp(— ^e) and the effect may have 
either sign depending on the sign of ^e- The most accurate recent bounds on 
the values of the chemical potentials based on a combined analysis of the angular 
spectrum of CMBR and BBN with the Deuterium fraction D/H = (3.0 ±0.4) ■ 10-^ 
are 

-0.01<ee<0.2, I^M.rl <2.6 (8) 

The third possible phenomenon is a non-equilibrium distortion of neu- 
trino spectrum. By the reasons specified above it has an especially strong effect on 



1 + 0.3488^ I 2 




(7) 
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Figure 4: Mass fraction of '^He as a function of the number of massless neutrino 
species. Different curves correspond to different values of the baryon-to-photon ratio 
r]io = lO^^ns/n^ = 2 (solid), 3.1 (dotted), 4 (dashed), and 5 (dashed- dotted). 



BBN for the case of electronic neutrinos. Normally neutrinos are well in thermal 
equilibrium which is not distorted by the cosmological expansion even if neutrino 
interactions is switched off at low temperatures, T < (1 — 2) MeV. The latter is true 
only for massless (or very light) particles. The spectrum may be distorted if there 
are some new massive particles (or some neutrinos are heavy) which can decay or 
annihilate into light neutrinos after their decoupling. The effect could be large but 
even in the standard model with the usual massless neutrinos a significant spectrum 
distortion exists IH' EID. Initially at large temperatures neutrinos were in strong 
thermal contact with electron-positron plasma, but when the temperature dropped 
down the contact became weaker and neutrino and electromagnetic components of 
the plasma became almost independent. Later e^e~ annihilation at T < me heats up 
the electromagnetic component of the plasma and the temperatures of neutrinos and 
electrons become different. Residual annihilation of hotter electrons, e^e~ vv, 
distorts neutrino spectrum. The effect was accurately calculated by the numerical 
solution of the integro-differential kinetic equations. The most precise calculations 
performed by two different groups and by different methods show a very 



good agreement. According to their results an excess of the energy densities for z/g 
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Figure 5: Deuterium-to-hydrogen by number as a function of the number of massless 
neutrino species. Notations are the same as in fig. 



and z/„ and z/,- are respectively 



hujpu^ = 0.94% and 6p./p. = 0.40% 



(9) 



However the impact of this effect on ^He is extremely small, at the level of lO""^. 
On the other hand, such correction may be in principle observed in high precision 
measurements of CMB anisotropics by future PLANCK satellite mission 



22, m 



if the canonical model can be tested with the accuracy of about 1% or better. A 
change in neutrino energy density with respect to the standard case would result in 
a shift of equilibrium epoch between matter and radiation, which is imprinted on 
the form of the angular spectrum of fluctuations of CMBR. The total energy density 
of relativistic matter in the standard model is given by 



rel 



1 + 0.68 — 



/1.401T, 



(10) 



where is the relative energy density of cosmic electromagnetic background radi- 
ation (CMBR) and is the photon temperature. The corrections discussed in this 
section and electromagnetic corrections of ref. HI) could be interpreted as a change 
of from 3 to 3.04. 



3 Effects of neutrino oscillations on BBN 



In the case of oscillations between active neutrinos, i.e. v^, i/^, and there would be 
no impact on BBN, if neutrinos were initially in thermal equilibrium with vanishing 
chemical potentials. Kinetic cqTiilibrium is established automatically at T > 2 
MeV due to a large weak interaction rate in comparison with the expansion rate. 
Vanishing of neutrino chemical potentials is the standard assumption but it is not 
necessarily true. If neutrinos are degenerate, oscillations between different flavors 
could have an impact on the production of light elements. The result depends 
not only on the parameters of the oscillations but also on the initial values of the 
lepton asymmetries in different sectors. More interesting physical effects originate 
in the case of oscillations between active and sterile {ug) neutrinos. Initially 
Ug should be absent in the primeval plasma. They can be produced only through 
the mixing with active ones but the mixing in matter is strongly suppressed at high 
temperatures. 

All three effects discussed in the previous section through which neutrinos 
could influence BBN, are operative in the case of oscillations between active and 
sterile neutrinos: 1) new relativistic species (namely Vg) are created, enlarging N,,] 
2) in the case of resonance oscillations a large lepton asymmetry may be generated 
in the sector of active neutrinos; 3) spectrum of active neutrinos, especially of i^e, 
can be distorted by oscillations. 

If there exist one or several types of sterile neutrinos, i.e. neutrinos which 
do not any interaction with other particles and manifest themselves only through 
non-diagonal mixing in the mass matrix, the eigenstates of the free Hamiltonian can 
be parametrized in the following way: 

Vl = Va COS 9yac + Sin 9yac, = -l^a Sm O^ac + l^s COS 9yac (H) 

where 9^ac is the mixing angle in vacuum. The effects of medium are described by 
the refraction index (or, what is practically the same, by effective potential). The 
latter can be calculated from the the Dirac equation averaged over thermal bath of 
cosmic plasma. The equation, roughly speaking, has the form: 

= ^ ij^dOaiJb) Oai^a (12) 

where Oa — 7a(l + 75)- The factor in the r.h.s. in front of ■0a, averaged over plasma, 
gives the effective potential, {Gp/ \/^){i^Oai^) — SatVeff- The space component of 
this operator vanishes due to isotropy of the plasma, while the time component is 



non-zero if the plasma is charge asymmetric. For the normally assumed magnitude 
of any cosmological charge asymmetry at the level 10^^ — 10^^° this term is sub- 
dominant and the main contribution to refraction index comes from the so-called 
non-local term. 

The operators ipa in the l.h.s. and r.h.s. of eq. (|T2|) may be at different 
space-time points because the weak interaction is not a local one but is mediated by 
an exchange of intermediate W or Z bosons. So the interaction is not the product 
of currents taken in the same space-time point. The contribution of non-locality is 
inversely proportional to m'^^ but still of the first order in the coupling constant. 
Taking all contributions together we obtain EJ): 

Veff = ±Cir]GFT^ + C^GlT^Ea-^ (13) 

where E is the neutrino energy, a = 1/137, Ci = 0.95, = 0.61, and C2'^ = 0.61. 
The numerical values of the coefficients Cj are calculated for thermally equilibrium 
bath and for temperatures below 100 MeV when muons were not present in the 
plasma. The coefficient rj is the charge asymmetry of the cosmological plasma, 
including all particle species interacting with z/q: 

= 2r/,^ + 7],^ + 7],^ +7],- 7]„J2 (for u,) , (14) 
r/(^) = 2r],^ + T],^ + T],^ - 7] J 2 (for v^) , (15) 

and T]'-'^^ for Uj. is obtained from eq. (|T5|) by the interchange /i ^ r. The indi- 
vidual charge asymmetries, rjx-, are defined as the ratio of the difference between 
particle-antiparticle number densities to the number density of photons: r]x = 
{nx — fix) /n^. The charge asymmetry terms comes with negative sign for neu- 
trinos and with positive sign for antineutrinos. 

At small T the charge asymmetry term may be larger than the non-local 
one but at smaller temperatures the vacuum term, dm^ cos29^ac/2E dominates for 
> 10~^ eV^. So for the "normal" value of the cosmological charge asymmetry 
the first term in eq. (|13|) is practically always sub dominant. 

An important difference between neutrino oscillations in stellar environ- 
ment and in cosmology is that in the former the loss of coherence is not important, 
neutrinos can be described by wave function and first order effects (in terms of the 
coupling constant) given by the effective potential are sufficient for description of all 
essential physics. Situation in cosmology is more complicated. Neutrino production 
and annihilation, as well as elastic scattering at non-zero angle could be quite strong 
at high temperatures, so the coherence is quickly destroyed and the density matrix 



formalism should be used il' 0' 113' HI) . The kinetic equation for density matrix 
in the cosmological background has the form 

P={j^- Hp^^ p = i[nm + Vess. p]+J driu, I, I) I^MAA+ - \ [p, ApA^}^ + 

j dril, v', I') [fvBp'B+ - i/, {p, BB+}^ (16) 

where p is the density matrix of the oscillating neutrinos, /; is the distribution 
function of other leptons in the plasma, dr is the phase space element of all particles 
participating in the reactions except for the neutrinos in question. It has essentially 
the same form as the usual collision term in kinetic equation for the number density, 
when oscillations are absent: 

If' = d^zdvy5\p,+py-pz)[\A{l + Y ^Z)\' 

n / 11(1 ± /)- - ^ + I' n / n(i ± /)] (i?) 

Y Z Z i+Y 

Here Y and Z are arbitrary, generally multi-particle states, Hy / is the product of 
phase space densities of particles forming the state F, and 

The signs '+' or '— ' in n(l i /) are chosen for bosons and fermions respectively. 

The first commutator term in the r.h.s. is first order in the interaction. 
It is the usual contribution from refraction index which does not break coherence. 
The last two terms are second order in the interaction and are related respectively 
to annihilation, uv //, and elastic scattering, vl v'V . The quantum statistics 
factors, (1 — /), and (/ — p) are neglected here. They can be easily reconstructed, 
see ref. In the interaction basis and in the case of active-sterile mixing the 

matrices A and B have only one non-zero entry in the upper left corner equal to 
the amplitude of annihilation or elastic scattering respectively; the upper means 
Hermitian conjugate. 

The exact form of the coherence breaking terms presented above is quite 
complicated. In many cases an accurate description can be obtained with a much 
simpler anzats: 

p = ...-{r,(p-pe,)} (19) 



where the multi-dots denote contributions from the neutrino mixing in medium de- 
scribed by the commutator term in eq. (0), Peq is the equilibrium value of the density 



matrix, i.e. the unit matrix multiplied by the equilibrium distribution function (^), 
and the matrix F that describes the interaction with the medium, is diagonal in the 
flavor basis; it is expressed through the reaction rates. If we take for the latter the 
total scattering rate, including both elastic scattering and annihilation, we obtain in 
the Boltzmann approximation ^ID: 



To = 2T,=g^-^GlT'p. (20) 

In general the coefficient ga{p) is a momentum-dependent function, but in the ap- 
proximation of neglecting [1 — /] factors in the collision integral it becomes a con- 
stant HO) equal respectively to g^^ ~ 4 and gi, ~ 2.9 0). In ref. 



more 



accurate values are presented that are found from the thermal averaging of the 
complete electro- weak rates (with factors [1 — /] included), which were calculated 



numerically from using the Standard Model code of ref. S). This gives g^^ ~ 3.56 
and g^^^^,^ ^2.5. 

Let us first consider non-resonance oscillations. The effective mixing angle 
in the medium is given by 

. asin2^^„e .^.x 
sin20,„ = — r— 21 

acos26'^„c - Veff 

where a = 5rn? /2E. In the pioneering papers 0) and in many subsequent ones 
the production rate of sterile neutrinos was estimated as: 

r, = (r,sin2 2^^sin2^) (22) 

where is the interaction rate of active neutrinos and after averaging over time one 
can substitute {sm^i^drn? t/2E)) = 1/2. Demanding that the energy density of z/^ 
produced by oscillations does not exceed the energy corresponding to AN^ neutrino 
species, one can obtain the following limit 

6m^ sin^ 9,ac < 6 ■ 10"^ (^i^) eV\ (23) 



3MeV/ 

where Ta is the decoupling temperature of active neutrinos, z/^. In ref. EIP the 
latter was taken 3 MeV for and 5 MeV for It corresponds to account only 
of the anmhilatiou rate. In ref. » the deeouphng ten,pe'atures were T, = 1.8 
eV and T^^.^ = 2.7 eV, corresponding to account of the total reaction rate. To 
resolve the problem one has to turn to the kinetic equations for the density matrix 
with the exact form of the coherence breaking terms (|TB|). Introducing real and 



imaginary parts for the non-diagonal elements of the density matrix according to 
Pas = P*sa = R + il, one finds: 

paaiPl) = -FI - / All \paa{pi)fl{p2) - Paa{p?i) fl{Pi)] - 



Knn [paa{Pl)paa{P2) " f l{Pz) fliP i)] , (24) 
PssiPl) = FI, (25) 



i?(j9i) = WI- (1/2) R{p 



rfr(z/2, I3, h)A^^^Paa{P2) ] , (26) 

/(pi) = -WR - (F/2) (p,, - p,a) - (1/2) /(pi) [I dr(/2, z/3, ^4)^2^,(^2) + 

dr(z/2, ^3, ^4)^LnPaa(P2) ] , (27) 

The integration is taken over the phase space according to: 

, 1 /■ d^P2 d^P2 d^p2 , .44. . , . 

'^'=2eJ (2.r)32E2 (2^)32^2 (2^)32^2 ^'"^^ (/>i + - - P4) (28) 
The amplitude of elastic scattering and annihilation with proper symmetrization 
factors can be taken from tables of ref. 

We will formally solve equations (p6| , |27|) to express the real and imaginary 
parts R and / of the non-diagonal components through the diagonal ones. The 
relevant equations can be written as 

R = WI-ToR/2, (29) 

/ = -WR+^{paa-Pss)-roI/2. (30) 

In the limit when the oscillation frequency 

u;,sc={F' + W'y^" (31) 

is much larger than the expansion rate, the solution is given by stationary point 
approximation, i.e. by the condition of vanishing the r.h.s.: 

FW 

^ ^ 2 {W^ + r2/4) ~ ^''^ ^^^^ 
FTo 

A{w^ + rl/A) (33) 

In the non-resonant case, when W ^ 0, usually the condition ^ is fulfilled 
and 

R ^ {sm29j2) {paa - Pss) (34) 

I ^ (sin 2^^ro/4l^) {paa - Pss) (35) 



where 6m is the mixing angle in matter and in the hmit of small mixing 

tan 2^^ ^ sin 2^^ ~ T77 = , ^oft/ — (36) 

W cos2d + {2EVeff/6m^) 

Now we can substitute the expression for I into eqs. (p^ . |25| ) and obtain a closed 
system of equations for the two unknown diagonal elements of density matrix which 
is easy to integrate numerically. In the case that the number density of sterile 
neutrinos is small and the active neutrinos are close to equilibrium, as often the case, 
we obtain the following equation that describes the production of sterile neutrinos 
by the oscillations 

p,, ^(sin2 2^^ro/4)/e, (37) 

An important conclusion of this derivation is that this production rate is 
by factor 2 smaller than the approximate estimates used in practically all earlier 
papers. An explanation of this extra factor 1/2 is that the time derivative of pss is 
proportional to imaginary part of the non-diagonal component of the density matrix 
and the latter is proportional to Fi = ro/2. 

The equations for the diagonal elements of density matrix can be solved 
analytically in the limit of small number density of sterile neutrinos, <^ Paa 
and for the case of kinetic equilibrium of active neutrinos, so that in Boltzmann 
approximation their distribution takes the form: 

Paa = C{x) exp(-y) (38) 

where dimensionless variables x = 1 MeVa(t) and y = pa{t) are introduced, with 
a{t) being the cosmological scale factor, normalized at high T as a = 1/T. Here 
C{x) can be understood as an effective chemical potential, C{x) = exp[,^(a;)], same 
for u and u. A justification for this approximation is a much larger rate of elastic 
scattering (that maintains the form ( pS] ) of Paa) with respect to annihilation rate 
that forces ^ down to zero (or C, = — iii the case of non-zero lepton asymmetry). 
This approach is similar to the standard calculations of cosmological freezing of 
species. Now we can integrate both both sides of eq. (0) over d^y so that the 
contribution of elastic scattering disappears and the following ordinary differential 
equation describing evolution of C{x) is obtained: 



dC _ h 
dx x'^ 



^ ^ + ^ 24(1 + 2gl + 2gl) J ^'"^ 



(39) 



The first term in the r.h.s. of this equation came from annihilation and the second 
one from the oscillations; the contribution of elastic scattering disappeared after 
integration over d^yi. 

We have assumed above that sin 29m <^ 1 and thus the term ~ (sin 29m)'^dC/dx 
has been neglected. It is a good approximation even for not very weak mixing. The 
constants ki are given by 

so that ke = 0.17 and kf^^r = 0.098. The charge asymmetry term in neutrino 
refraction index was also neglected. It is reasonable in the non-resonant case if the 
asymmetry has a normal value around 10^^ — 10^^°. 



Eq. (|39D can be solved analytically if|5| = |l — C|^l and after some 
algebra we obtain the following result for the increase of the effective number of 
neutrino species induced by mixing of active neutrinos with sterile ones: 



where 



1 sin' 26,,, GUI + 9l + 9'r) 9*iT;^odJ .... 
^"97r2 ^ Hx^ 10.75 ^ ' 



2.34 ■ 10-^ H /9 - 0-65 • 10"^ . . 

~ 5m2 cos 29 ^^'^ ~ 5w? cos 26 ' ^ ' 

Substituting numerical values of the parameters we obtain: 

(5<../eV2) sin^2e:;= = 3.16 ■ 10-^(^?.(T-„,)/10.75)3(AiV.)2 (43) 

('^<../eV2) sin' 2ec" = 1-74 ■ 10"^(^?.(T- )/10.75)3(AiV.)2 (44) 



Here another factor g^, came from the Hubble parameter. 

These results can be compared with other calculations. They are ap- 
proximately 2 orders of magnitude stronger than those presented in refs. 0), 
where too high freezing temperature for weak interaction rates was assumed and 
the limit was obtained: (5m^sin^2^ < 6 ■ 10~^AA^,^. In ref. 0) the limit was 



(5m^sin^2^ < 3.6 ■ 10 ^AN'^. (All these are given for mixing with Ue-) On the other 
hand, the limits obtained in ref. are approximately 6 time stronger than those 
found above (|^P^ ). They are: sin^2^ < 5 ■ 10~^AA''^ for mixing with z/g and 
(5m^sin^26' < 3 ■ lO'^AA'^^ for mixing with u^^r- The difference by factor 6 between 
these results and eq. ( |43|J4^ ) can be understood in part by the factor 2 difference 



in the interaction rate, according to eq. {^7\), which gives factor 4 difference in the 
limits. The remaining difference by roughly factor 1.5 could possibly be prescribed 



to different ways of solution of kinetic equations or to the fact that the increase in 
the number density of sterile neutrinos is accompanied by the equal decrease in the 
number density of active neutrinos if the production of the latter by inverse anni- 
hilation is not efficient. This phenomenon is missed in kinetic equations with the 
simplified form of the coherence breaking term (|T^, which are mostly used in the 
literature, while equations (P^-P7|) automatically take that into account. However, 
for sufficiently large mass difference, 5rn?, the effective temperature of Vs production 
is: 

= (12 - 15) (3/^)1/=^ (5mVeV2)V6 MeV (45) 

Here the first number is for mixing of Vs with and the second one is for mixing 
with z/^ or Vt-. It is larger than the temperature of the annihilation freezing, so 
the active neutrino states are quickly re-populated and the mentioned above affect 
could be significant only for a small mass difference. Much weaker bounds obtained 
in ref. 0) resulted from an error in the coherence breaking terms in kinetic equations 
for the non-diagonal matrix elements of the density matrix. 

The result (|i3| ) includes only a rise of the total cosmological energy density 
due to production of sterile neutrinos. In the case of mixing with z/g the distortion 
of spectrum of the latter is also important for BBN. This issue is addressed in 
refs. H, i3 P. 



If the mass difference between and z/^ is negative then the MSW-resonance 
transition might take place in cosmological background. The resonance may also 
exist with an arbitrary sign of 5m? if the initial value of the asymmetry is suffi- 
ciently large but we will not consider this case. A very interesting phenomenon 



was discovered in ref. 0). It was found that initially very small lepton asymme- 
try (~ 10"^) can rise almost to unity because of oscillations. The effect is based 
on the dependence of the refraction index on the cosmological charge asymmetry 



(the first term in eq. ([T3|) ) . Since this term comes with different signs to neutrinos 
and antineutrinos, the resonance conditions for the same value of momenta could 
be reached e.g. earlier for v than for v. This would give rise to a more efficient 
transformation of neutrinos in comparison with antineutrinos and the asymmetry 
could start rising. This creates a positive feed-back effect and exponential rise of 
the asymmetry at initial stage. In fact, initially the asymmetry very strongly drops 
down and only after some rather short period the explosive rise begins. The first 
results were obtained by numerical solution of kinetic equations for the density ma- 
trix, with some reasonable approximations made. In subsequent works more precise 
numerical methods have been developed. The present-day state of art and the list 



T/MeV 



Figure 6: Evolution of L^^^ = 2L^^ + rj for v^. — > z/^ oscillations with sin^ 29o = 
10~^ and, from left to right, dm? /eV"^ = —0.25, —0.5, —1.0, —2.0, —4.0 obtained from 
numerically solution of the quantum kinetic equations. The initial L^^ = is taken 
and rj = 5 X 10~^° is assumed. The low temperature evolution is weakly dependent 
on these values. 



of relevant references can be found e.g. in the recent papers 0' © clP. Numerical 
solution presents a serious challenge because of quickly oscillating functions under 
the sign of the integral over neutrino momentum which enters the charge asymmetry 
term. On the other hand, the presence of fast and slow variables permits to separate 
them and to find an approximate analytical solution ED . Numerical and analytical 
results are in a very good agreement in the coinciding range of parameters. The 
evolution of z/g-asymmetry according to calculations of ref. 0) for sin^ 26 = 10~^ 
and several values of mass difference is presented in fig. ^. 

Immediately after the large rise of lepton asymmetry due to oscillations 
between active and sterile neutrinos was discovered, it was argued that the 
resulting asymmetry is large but chaotic, i.e. its sign is essentially unpredictable. 
The sign of the asymmetry is very sensitive to the oscillation parameters and the 
input of numerical calculations. As a result of this feature the sign of lepton asym- 
metry might be different in different causally non-connected domains E). This 
could have interesting implications and, in particular, would lead to inhomogeneous 
nucleosynthesis. The existence of chaoticity was confirmed in several subsequent 



papers but again in the frameworks of simphfied thermally averaged equations. 
On the other hand, the analysis of possible chaoticity performed in ref. ED on the 
basis of numerical solution of kinetic equations with a full momentum dependence 
shows a different picture. Most of parameter space is not chaotic, while in the re- 
gion where chaoticity is observed numerical calculations are not reliable. Analytical 
calculations of ref. E3) also do not show any chaoticity. At the present time the 
problem remains unresolved. 

A very interesting phenomenon of spatial fluctuations of lepton asymme- 
try was found in ref. Neutrino oscillations in the presence of initially small 
baryonic inhomogeneities could give rise to domains with different signs of lepton 
asymmetry. This effect is different from chaotic amplification of asymmetry dis- 
cussed above. As we have already mentioned the initial asymmetry first drops down 
to an exponentially small value and after that started to rise, also exponentially, 
with a larger integrated exponent. Since the value of the asymmetry in the mini- 
mum could be extremely small, it is sensitive to small perturbations and the final 
sign of asymmetry could be determined by them. This goes not in a trivial way as 
e.g. spatial fluctuations of the sign in the minimum but somewhat more tricky and 
related to diffusion term in kinetic equations. As is argued in the original paper ED 
(see also ref. , small spatial fluctuations of the cosmological baryon asymmetry, 
though they do not create change of sign in the minimum, would induce formation of 
domain with super-horizon sizes (at the moment of their creation) with large lepton 
asymmetry of different signs. 

Let us consider the equation, modeling evolution of the lepton asymmetry 
in the presence of small spatial inhomogeneities, used in ref. 0). Following nota- 
tions of this paper let us denote the asymmetry of active neutrinos of flavor a as 
Ly^. The combination that enters the refraction index of Ua is 2L,y^ + L, where L 
consists of the contributions of other neutrino species, baryons, and electrons. Due 
to electric charge neutrality the last two are not independent. The preexisting fluc- 
tuations in asymmetry of neutrinos would be erased due to a large neutrino mean 
free path in cosmic plasma after neutrino decoupling. So the background asymmetry 
can be written as L = L + 6B{x), where the first term is homogeneous and the last 
one could be inhomogeneous and related to the fluctuations in the baryon number. 
Now we will neglect the universe expansion (it is essential for quantitative estimates 
but not for a qualitative result) and will solve the evolution equations in somewhat 
different way, permitting more accurate evaluation of the effects of diffusion. The 



evolution of the neutrino asymmetry is described by the equation: 



(f , t) = a{t) [2 L,^ {x, t) + L + SB{x)] + D{t)V''L,^ {x, t) (46) 

where D(t) is the diffusion coefficient, and the function a{t) is initiaUy negative and 
generate an exponential decrease of the asymmetry, but at some critical time tc it 
changes sign and this creates a huge rise of the asymmetry. It is essential that, while 
a{t) is negative, the asymmetry drops down to a very small value. In a more accurate 
formulation a would also depend on the asymmetry itself, but in what follows we 
are interested in rather small values of the asymmetry, where the non-linear effects 
are not important. 

The solution to this equation can be found by the Fourier transform and 
we obtain: 

{x, t) = L f dt' ait') C + (47) 

+ t dt'a{t')e^^>'"<'"^ I d'k(^^^5B{k,t,^)e-''"C^^"^^n 

Here "hut" indicates the Fourier transform of the corresponding function. The first 
integral in this expression can be explicitly taken because the integration measure 
dt'a{t') is exactly the differential of the exponential. The integration of this term 
gives: 



(1/2)Z 



exp (2 f dt2a{t2)] - 1 



(48) 



So we obtained a rising term (after some initial decrease) plus a constant initial 
value of L. 

Moreover, the second term can be further integrated because the initial 
value L{k,tin) is supposed to be homogeneous and so its Fourier transform is just 
delta-function, S^{k). The integral gives 

exp (2 dt2a{t2)) (49) 

V *' tin ^ 

where L^^"^^ means the initial value, i.e. taken t — tin- So if we forget about 
the constant term Z/2, we would have {L^J^^ + L/2) multiphed by the rising expo- 
nent. One would get exactly this expression if one solves the equation for L^,^ in 
homogeneous case. 



Let us now consider the last term which, according to the arguments of 
ref. could change the sign of the rising asymmetry, i.e. this last term could 

become the dominant one. To evaluate the integral let us substitute: 

SB{k,ti^) = f d^xie''''^' 6B{xi) (50) 



where 6B{xi) is the initial value of the inhomogeneous term. Now we can make 
integration over d^k. We have the integral of the type 

J d^kexpl-S'^k'^ + ik{x-xi)] (51) 

the scalar product of vectors k and r = x — xi is equal to k{x — xi) = kr cos 9, and 

52(t) = /■* dhD{h) (52) 

Integration over angles in d^k = 2TTk'^dkd{cos6) is trivial, it gives sinkr/kr. The 
remaining integration can be done as follows: 

y dkk sin kr exp[—S'^k^] = [d/dr) j dk cos kr exp[—S'^k'^] (53) 

and the last integral can be taken if we expand the range of integration from minus 
to plus infinity. Introducing new variable Xi = x — S{ti)pwe finally obtain 

J dtia{h) e ^''"^''^ j d'p 6B{x - S{ti)p) e'^' (54) 

This is the contribution the lepton asymmetry Lj^^ generated by the (small) baryonic 
inhomogeneities. Its asymptotic rise at large t is similar to the rise of other terms 
but its exponential decrease at intermediate stage could be considerably milder and, 
as a result, this term could become dominant with the sign determined by the sign 
of the fluctuations in the baryon asymmetry. We can check this on a simple example 
assuming that the function a{t) has the form a{t) = ai{t — t^) and that fluctuations 
of the asymmetry are described by one harmonic mode: 6B{x) = escoskox. This 
form of SB could be inserted either into eq. ( p^ or into initial eq. ( ^8D and we find 
for the oscillating part of the asymmetry (up to a constant coefficient): 

SL{x) = eBCOsfcoxe'^^(*-*=)'-^W'=o [l*"*^ rft,t,e-"^*?+^(*^)'=o + 

(55) 
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Both terms rise as exp[ai(t — tc)^], i-e. in the same way as the other homogeneous 
terms (we assume that S{t) is finite at large t and not too large). The first term is 



exponentially suppressed as exp[— ai(tc — ^m)] also at the same level as the homo- 
geneous terms. The second term, which vanishes in homogeneous case ( fco = or 
5* = 0) is not exponentially suppressed. In the limit of large ai the integral can be 
evaluated as ~ >S'(0) kQ/a\. It is small but not exponentially small. Thus, it is easy 
to imagine the situation when the last term dominates and resonance enhancement 
of lepton asjTiimetry in the background of small fluctuations of baryon asymmetry 
could create domains with large and different lepton asymmetry. The effect is very 
interesting and deserves more consideration. 

4 Heavy sterile neutrinos (with 10-100 MeV mass) 

A hypothesis that such neutrinos may exist originated from the observation of the 
KARMEN anomaly in the time distribution of the charged and neutral current 
events induced by neutrinos from 7r+ and /i"*" decays at rest A suggested 

explanation of this anomaly was the production of a new neutral particle in pion 
decay 

7r+ ^ /i+ + x° , (56) 

with the mass 33.9 MeV, barely permitted by the phase space, so that this parti- 
cle moves with non-relativistic velocity. Its subsequent neutrino-producing decays 
could be the source of the delayed neutrinos observed in the experiment. Among 
possible candidates on the role of x'^-particle was, in particular, a 33.9-MeV sterile 
neutrino 0) . According to ref. EH) , cosmology and astrophysics practically exclude 



the interpretation of the KARMEN anomaly by a 33.9 MeV neutrino mixed with v-r. 
It agrees with a later statement by the KARMEN collaboration made at Neutrino 
2000 ^3) that the anomaly was not observed in upgraded detector KARMEN 2, 



but the question still remains which area in the mass-mixing-plane for heavy sterile 
neutrinos can be excluded. This issue was addressed recently by NOMAD collab- 
oration in direct experiment ^ID and in ref. 0) by considerations of big bang 
nucleosynthesis and the the duration of the supernova (SN) 1987A neutrino burst. 

The mixing of heavy sterile neutrinos with the active ones couples the 
former to boson and induce the decay 

V2-^vi + l + I, (57) 

where i is any lepton with the mass smaller than the mass m2 of the heavy neutrino. 
If m2 < 2m^ the decay into and fr is kinematically forbidden. If Vs is mixed 



either with z/^ or z/,-, the hfe-time is expressed through the mixing angle as: 



~ (M,/10 MeV)5 sin22^ ' ^ ' 

For the mixing with z/g the numerator is 0.7 sec; the difference is due to the charged- 
current interactions. 

There are several effects operating in different directions, by which a heavy 
unstable sterile neutrino could influence big-bang nucleosynthesis. First, their con- 
tribution to the total energy density would speed up the expansion and enlarge the 
frozen neutron-to-proton ratio. Less direct but stronger influence could be exerted 
through the decay products, z/g, z/^, and Ur, and and through the change of the 
temperature evolution, T^/T^. The impact of z/^ and on BBN is rather straight- 
forward: their energy density increases with respect to the standard case and this 
also results in an increase of r„. This effect can be described by the increased num- 
ber of effective neutrino species Ny during BBN. The increase of the energy density 
of z^e, due to decay of z/^ into z/g, has an opposite effect on r„. Though a larger energy 
density results in faster cooling, the increased number of z/g would preserve thermal 
equilibrium between neutrons and protons for a longer time and correspondingly 
the frozen n/p-ratio would become smaller. The second effect is stronger, so the 
net result is a smaller n/p-ratio. There is, however, another effect of a distortion 
of the equilibrium energy spectrum of z/g due to produced from the decays of z/^. 
If the spectrum is distorted at the high-energy tail, as is the case, then creation of 
protons in the reaction + z/g — > p + e~ would be less efficient than neutron creation 
in the reaction z/g + p — ^ n + e"^. We found that this effect is quite significant. Last 
but no the least, the decays of Vg into the e^e~-channel will inject more energy into 
the electromagnetic part of the primeval plasma and this will diminish the relative 
contribution of the energy density of light neutrinos and diminish r„. 

In refs. S' the distribution functions of neutrinos were calculated 
from kinetic equations in Boltzmann approximation and in a large part of parame- 
ter space they significantly deviate from equilibrium. The distributions of electrons 
and positrons were assumed to be very close to equilibrium because of their very 
fast thermalization due to interaction with the photon bath. However, the evolu- 
tion of the photon temperature, due to decay and annihilation of the massive Vg 
was different from the standard one, Ty ~ l/a(t), by an extra factor (1 + A) > 1 
and this was explicitly calculated from the energy balance condition. At sufficiently 
high temperatures, T > Tw ~ 2 MeV, light neutrinos and electrons/positrons were 
in strong contact, so that the neutrino distributions were also very close to the 



equilibrium ones. If Ug disappeared sufficiently early, while thermal equilibrium be- 
tween and neutrinos remained, then z/^ would not have any observable effect on 
primordial abundances, because only the contribution of neutrino energy density 
relative to the energy density of and 7 is essential for nucleosynthesis. Hence a 
very short-lived has a negligible impact on primordial abundances, while with an 
increasing lifetime the effect becomes stronger. Indeed at T < Tw the exchange of 
energy between neutrinos and electrons becomes very weak and the energy injected 
into the neutrino component is not immediately redistributed between all the parti- 
cles. The branching ratio of the decay of Ug into e'^e~ is approximately 1/9, so that 
the neutrino component is heated much more than the electromagnetic one. As we 
mentioned above, this leads to a faster cooling and to a larger ra/p-ratio. 

If the equilibrium number density of sterile neutrinos is reached, it would be 
maintained until Tf ^ 4(sin2^)~^/^ MeV. This result does not depend on the heavy 
neutrino mass because they annihilate with massless active ones, 1^2 + ^0 cd^- The 
heavy neutrinos would be relativistic at decoupling and their number density would 
not be Boltzmann suppressed if, say, Tf > Ms/2. This gives 

sin^ 2d{5m^ /MeY^f^ < 500 . (59) 

If this condition is not fulffiled the impact of on BBN would be strongly dimin- 
ished. On the other hand, for a sufficiently large mass and non-negligible mixing, 
the z/2 hfetime given by Eq. ( |5^ ) would be quite short, so that they would all decay 
prior to the BBN epoch. (To be more exact, their number density would not be 
frozen, but follow the equilibrium form oc e~"^^^'^~'.) 

Another possible effect that could diminish the impact of heavy neutrinos 
on BBN is entropy dilution. If z/2 were decoupled while being relativistic, their 
number density would not be suppressed relative to light active neutrinos. However, 
if the decoupling temperature is higher than, say, 50 MeV pions and muons were still 
abundant in the cosmic plasma and their subsequent annihilation would diminish 
the relative number density of heavy neutrinos. If the decoupling temperature is 
below the QCD phase transition the dilution factor is at most 17.25/10.75 = 1.6. 
Above the QCD phase transition the number of degrees of freedom in the cosmic 
plasma is much larger and the dilution factor is approximately 5.5. However, these 
effects are essential for very weak mixing, for example the decoupling temperature 
would exceed 200 MeV if sin^26' < 8 x 10~^. For such a small mixing the life-time 
of the heavy 1^2 would exceed the nucleosynthesis time and they would be dangerous 
for BBN even if their number density is 5 times diluted. 

Sterile neutrinos would never be abundant in the universe if Tg/H < 1. In 



r 



5.8- 10"^ sec 



■^.^2 ^m^(m2 - ) 



(61) 



fact we can impose a stronger condition demanding that the energy density of heavy 
neutrinos should be smaller than the energy density of one light neutrino species at 
BBN (T ~ 1 MeV). Taking into account a possible entropy dilution by factor 5 we 
obtain the bound: 

((JmVMeV^) sin^ 26 < 2.3 x 10"^ . (60) 

Parameters satisfying this conditions cannot be excluded by BBN. 

If Vs (or to be more precise 1^2) mass is larger than 135 MeV, the dominant 
decay mode becomes z/2 — * vr*^ + z/^. The life-time with respect to this decay can be 
found from the calculations 0' 0) of the decay rate 7r° — > vi) and is equal to: 

Gpms{ml -ml)f^ sin^ O]'^ 
167r 

where = 1712 is the mass of the sterile neutrino (we always assume that the mixing 
angle is small so that 1^2 ~ i^s), = 135 MeV is the 7r°-mass and /^r = 131 MeV 
is the coupling constant for the decay 7r+ fi + u^. The approximate estimates 
of ref. permit one to conclude that for the life-time of 1^2 smaller than 0.1 
sec, and corresponding cosmological temperature higher than 3 MeV, the decay 
products would quickly thermalize and their impact on BBN would be small. For a 
life-time larger than 0.1 sec, and T < 3 MeV, one may assume that thermalization 
of neutrinos is negligible and approximately evaluate their impact on BBN. If 
is mixed with z/^ or then electronic neutrinos are not produced in the decay 
Ug Ti^Va and only the contribution of the decay products into the total energy 
density is essential. As we have already mentioned, non-equilibrium z/e produced by 
the decay would directly change the frozen n/p-ratio. This case is more complicated 
and demands a more refined treatment. 

The 7r° produced in the decay z/^ — > z/^ + vr'' immediately decays into two 
photons and they heat up the electromagnetic part of the plasma, while neutrinos by 
assumption are decoupled from it. We estimate the fraction of energy delivered into 
the electromagnetic and neutrino components of the cosmic plasma in the instant 
decay approximation. Let = ng/riQ be the ratio of the number densities of the 
heavy neutrinos with respect to the equilibrium light ones, no = 0.09T^. The total 
energy of photons and e+e^-pairs including the photons produced by the decay is 

T^ + r,no-^ 1 + ^ 
2 30 2 \ mi 



Pem = — TTTT^ + ^^^0^ I 1 + ri ' (62) 



while the energy density of neutrinos is 



Pu = + rsUo— 1 . 63 

4 30 2 I / ^ ' 



The effective number of neutrino species at BBN can be defined as 

Nlf"^ = f ^ . (64) 

' Pem 

Because of the stronger heating of the electromagnetic component of the plasma by 
the decay products, the relative role of neutrinos diminishes and Nl^^^") becomes 
considerably smaller than 3. li Vg are decoupled while relativistic their fractional 
number at the moment of decoupling would be = 4 (two spin states and antipar- 
ticles are included). Possible entropy dilution could diminish it to slightly below 1. 
Assuming that the decoupling temperature of weak integrations is Tw — 3 MeV we 
find that N^'^^f'^ = 0.6 for m, = 150 MeV and N^^^f'^ = 1.3 for = 200 MeV, if the 
frozen number density of Vg is not diluted by the later entropy release and remains 
equal to 4. If it was diluted down to 1, then the numbers would respectively change 
to A^^^-^-^) = 1.15 for rus = 150 MeV and N^^^^^ = 1.7 for = 200 McV, instead of 
the standard Njf^^^ = 3. Thus a very heavy would result in under-production of 
'^He. There could, however, be some other effects acting in the opposite direction. 

Since z/g decouples from electrons/positrons at smaller temperature than 
z/^ and Ur, the former may have enough time to thermalize. In this case the tem- 
peratures of i>e and photons would be the same (before e^e~-annihilation) and the 
results obtained above would be directly applicable. However, if thermalization be- 
tween Ue and was not efficient, then the temperature of electronic neutrinos at 
BBN would be smaller than in the standard model. The deficit of Ug would produce 
an opposite effect, namely enlarging the production of primordial ^ife, because it 
results in an increase of the n/p-freezing temperature. This effect significantly dom- 
inates the decrease of N^^f^^ discussed above. Moreover even in the case of the decay 
i/2 — > 7r° -|- Uf^^T, when i/g are not directly created through the decay, the spectrum 
of the latter may be distorted at the high energy tail by the interactions with non- 
equilibrium Ur and z/^ produced by the decay. This would result in a further increase 
of ^iJe-production. In the case of direct production of non-equilibrium through 
the decay z/2 — > 7r° -|- Ue their impact on n/p ratio would be even much stronger. 

To summarize, there are several different effects from the decay of into 
7r° and u on BBN. Depending upon the decay life-time and the channel these effects 
may operate in opposite directions. If the life-time of z/2 is larger than 0.1 sec 
but smaller than 0.2 sec, so that e"^ and Ug establish equilibrium the production 
of '^He is considerably diminished so that this life-time interval would be mostly 
forbidden. For life-times larger than 0.2 sec the dominant effect is a decrease of 
the energy density of Ug and this results in a strong increase of the mass fraction 
of ^He. Thus large hfe-times should also be forbidden. Of course there is a small 



part of the parameter space where both effects cancel each other and this interval 
of mass/mixing is allowed. It is, however, difficult to establish its precise position 
with the approximate arguments used in ref. . 

Thus, in the case of Us ^ v^^t mixing for nig > 140 MeV we can ex- 
clude the life-times of roughly larger than 0.1 sec, except for a small region 
near 0.2 sec where two opposite effects cancel and the BBN results remain undis- 
turbed despite the presence of sterile neutrinos. Translating these results into mix- 
ing angle according to eq. (^), we conclude that mixing angles sin^ ^ < 5.8 ■ 
10~^m^/ms/((ms/m^)^ — 1) are excluded by BBN. Combining this result with eq. 
(|60|) we obtain the exclusion region for > lAOMeV: 

5.1 . 10-8 < sin2 < 5.8 ■ 10-« . (65) 

mj nis [ms/mT^Y — 1 

In the case of z/^ ^ v^. mixing for rris > lAOMeV the limits are possibly 
stronger, but it is more difficult to obtain reliable estimates because of a strong 
influence of non-equihbrium i/g, produced by the decay, on neutron-proton reactions. 

The constraints on the mass/mixing of Ug from neutrino observation of 



SN 1987A are analyzed in some detail in ref. and based on the upper limit of 
the energy loss into a new invisible channel because the latter would shorten the 
neutrino burst from this supernova below the observed duration. 

The results are summarized in fig. 0. The region between the two horizontal 
lines running up to 100 MeV are excluded by the duration of the neutrino burst from 
SN 1987A. BBN excludes the area below the two upper dashed lines if the heavy 
neutrinos were abundant in the early universe. These two upper dashed lines both 
correspond to the conservative limit of one extra light neutrino species permitted 
by the primordial ^He-abundance. The higher of the two is for mixing with z/^^^ 
and the slightly lower curve is for mixing with Ue- In the region below the lowest 
dashed curve the heavy neutrinos are not efficiently produced in the early universe 
and their impact on BBN is weak. For comparison we have also presented here the 
region excluded by NOMAD Collaboration for the case of z/^ ^ mixing. A 
more accurate consideration would probably permit to expand the excluded region 
both in the horizontal and vertical directions. 



5 Acknowledgements 



I am greateful to P. Di Bari for the discussion of the formation of leptonic domains. 




00 







1 
10 











Excluded by BBN 
not excluded 



eV 



00 



Figure 7: Summary of the exclusion regions in the (siii^ 9-Ms)-plane. 



References 



1. A.T. Lee et al, |astro-ph/ 0104459 



C.B. Netterfield et al., |astro-ph/ 01 04460 



N.W. Halverson et al., ^stro-ph/ 0104489 



2. K.A. Olive, G. Steigman, and T.P. Walker, Phys.Rept. 333, 389 (2000). 

3. N. Hata, R.J. Scherrer, G. Steigman, D. Thomas, and T.P. Walker, Ap.J. 458, 
637 (1996). 



4. K.A. Olive, E. Skillman, and Gary Steigman, |astro-ph/9611166 . 

5. A. Peimbert, M. Peimbert, and V. Luridiana, astro-ph/010718S . 

6. R. Gruenwald, G. Steigman, and S.M. Viegas, |astro-ph / 1 090 7T 



7. T.X. Tliuan and Yu.I. Izotov, |astro-ph/ 0003234 



8. J.M. O'Meara, D. Tytler, D. Kirkman, et al, |astro-ph/0011179 



9. J.M. O'Meara, D. Tytler, D. Kirkman, et al, Astrophys.J. 552 718 (2001). 
10. S. D'Odorico, M. Dessauges-Zavadsky, and P. Molaro, |astro-ph/0 102162". 



11. M. Pettini and D.V. Bowen, |astro-ph/ 01 04474] ; Astrophysical Journal, scheduled 
for Vol. 560 (2001). 



12. P. Molaro, P. Bonifacio, M. Centurion, and G. Vladilo, |astro-ph/ 9908060 



13. S. A. Levshakov, M. Dessauges-Zavadsky, S. D'Odorico, and P. Molaro |astro- 
pli/0l0552-g|. 



14. J.K. Webb, R.F. Carswell, K. M. Lanzetta, et al. |astro-pli/97niUgU 



15. A.D. Dolgov and B.E.J. Pagel, New Astron. 4, 223 (1999). 



16. D. Kirkman, D. Tytler, J.M. O'Meara, et al. |astro-ph/01 03305 



17. S.H. Hansen, G. Mangano, A. Melchiorri, G. Miele and O. Pisanti, |astro- 



pli/0105385 



18. A.D. Dolgov and M. Fukugita, JETP Lett. 56 (1992) 123 [Pisma Zh. Eksp. 
Teor. Fiz. 56 (1992) 129]; 

A.D. Dolgov and M. Fukugita, Phys. Rev. D 46 (1992) 5378. 

19. S. Dodelson and M. S. Turner, Phys. Rev. D 46 (1992) 3372. 

20. A.D. Dolgov, S.H. Hansen and D.V. Semikoz, Nucl. Phys. B 503 (1997) 426 
||hep-ph/ 970331511; 



A.D. Dolgov, S.H. Hansen and D.V. Semikoz, Nucl. Phys. B 543 (1999) 269 
hep-ph/9805467| . 



21. S. Esposito, G. Miele, S. Pastor, M. Peloso and O. Pisanti, Nucl. Phys. B 590 
(2000) 539 ||astro-ph/000"5573[ . 

22. N.Y. Gnedin and O.Y. Gnedin, Astrophys. J. 509 (1998) 11 ||astro-ph/97r2T9S 
v3]. 

23. R.E. Lopez, S. Dodelson, A. Heckler and M.S. Turner, Phys. Rev. Lett. 82 
(1999) 3952 [|astro-ph/9803095| . 



24. D. Ndtzold and G. Raffeh, Nucl. Phys. B 307 (1988) 924. 

25. A.D. Dolgov, Yad.Fiz. 33 (1981) 1309; Enghsh translation: Sov. J. Nucl. Phys. 
33 (1981) 700. 

26. L. Stodolsky, Phys. Rev. D 36 (1987) 2273. 



27. G. Sigl and G. Raffelt, Nucl. Phys. B 406 (1993) 423. 

28. B.H.J. McKellar and M.J. Thomson, Phys. Rev. D 49 (1994) 2710. 

29. R.A. Harris and L. Stodolsky, Phys. Lett. B 116 (1982) 464. 

30. N.F. Bell, R.R. Volkas and Y.Y.Y. Wong, Phys. Rev. D 59 (1999) 113001 |Ei^ 



ph/9809363 |. 



31. K. Enqvist, K. Kainulainen and M. Thomson, Nucl. Phys. B 373 (1992) 498. 

32. A.D. Dolgov, S.H. Hansen, S. Pastor and D.V. Semikoz, Astropart. Phys. 14 



(2000) 79 [|liep-ph/99 10444 

33. R. Barbieri and A.D. Dolgov, Phys. Lett. B 237 (1990) 440. 
R. Barbieri and A.D. Dolgov, Nucl. Phys. B 349 (1991) 743. 

34. K. Kainulainen, Phys. Lett. B 244 (1990) 191. 

35. A.D. Dolgov, Phys. Lett. B 506 (2001) 7 ||hep-ph/0006T0^ . 

36. D.P. Kirilova and M.V. Chizhov, Phys. Lett. B 393 (1997) 375 |Ei^ 



ph/9608270 |. 



37. D.P. Kirilova and M.V. Chizhov, Phys. Rev. D 58 (1998) 073004 ^epj 
ph/970728^ 



38. D.P. Kirilova and M.V. Chizhov, Nucl. Phys. B 591 (2000) 457 jhe^ 



ph/9909408 



39. R. Foot, M.J. Thomson and R.R. Volkas, Phys. Rev. D 53 (1996) 5349 jhep^ 
ph/9509327 |. 



40. R. Foot and R.R. Volkas, Phys. Rev. D 61 (2000) 043507 ||hep-ph/990433^ . 

41. R. Foot, Phys. Rev. D 61 (2000) 023516 ||hep-ph/99063TT|] . 

42. P. Di Bari, R. Foot, |hep-ph/0103T92 . 

43. A.D. Dolgov, Nucl.Phys. B610, 411 (2001) [|hep-ph/0102125|| . 

44. P. Di Bari and R. Foot, Phys. Rev. D 63 (2001) 043008 ||hep-ph/0008"258[ . 

45. X. Shi, Phys. Rev. D 54 (1996) 2753 [ |astro-ph/960"2T35| . 



46. X. Shi and G.M. Fuller, Phys. Rev. Lett. 83 (1999) 3120 ||astro-ph/ 9904041 



47. K. Enqvist, K. Kainulainen and A. Sorri, Phys. Lett. B 464 (1999) 199 jEep^ 
ph/9906452| . 



A. Sorri, Phys. Lett. B 477 (2000) 201 |[hep-ph/99inUB |. 

P.N. Braad, S. Hannestad, |hep-ph/0012T94| . R. Buras, |hep-ph/ 0002086 . 

48. P. Di Bari and R. Foot, Phys. Rev. D 61 (2000) 105012 ||hep-ph/99im3 



49. P. Di Bari, Phys. Lett. B 482 (2000) 150 ||hep-ph/9911214 



50. K. Enqvist, K. Kainulainen and A. Sorri, JHEP 0104 (2001) 012 jhep^ 
ph/001229T 



51. B. Armbruster et al. [KARMEN Collaboration], Phys. Lett. B 348 (1995) 19; 
K. Eitel et al. [KARMEN Collaboration], Nucl. Phys. Proc. Suppl. 77 (1999) 
212 [ |hep-ex/980900l . 

52. V. Barger, R.J. Phillips and S. Sarkar, Phys. Lett. B 352 (1995) 365 [Erratum- 
ibid. B 356 (1995) 617] ||hep-ph/9503295|] . 

53. A.D. Dolgov, S.H. Hansen, G.G. Raffelt and D.V. Semikoz, Nucl. Phys. B 580 
(2000) 331 ||hep-ph/ 00022231 . 

54. K. Eitel [KARMEN Collaboration], Nucl. Phys. Proc. Suppl. 91 (2000) 191 
|hep-ex/0008002|] . 



55. P. Astier et al. [NOMAD Collaboration], Phys. Lett. B 506 (2001) 27 jh^ 
ex/0101041 



56. A.D. Dolgov, S.H. Hansen, G.G. Raffelt and D.V. Semikoz, Nucl. Phys. B 590 
(2000) 562 ||hep-ph/ 00081381 . 

57. R. Fishbach et al, Phys. Rev. D 16 (1977) 2377. 

58. T. Kalogeropoulos, J. Schechter and J.W.F. Valle, Phys. Lett. B 86 (1979) 72. 



